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Abstract 

We review the measurements of the top quark mass by the CDF and D0 collaborations 
using Run I data in excess of 100 pb“^. The D0 collaboration [P has recently updated 
its measurement of the top quark mass in the lepton + jets channel. The world average of 
the top quark mass from the CDF Q and D0 measurements in the lepton + jets channel 
now stands at 175.6 ± 5.5 GeV/c^. 


1 Introduction 


The top quark is of fundamental importance in understanding the standard model of particle 
interactions. The measurement of the top quark mass and the mass of the W boson can be 
used together to constrain the mass of the Higgs particle, an as yet unobserved particle which 
is responsible for the generation of masses in the standard model. Failure to obtain Higgs mass 
constraints within reasonable limits can become the hrst indication of the incompleteness of the 
standard model. This paper reports on the latest measurements of the top quark mass from 
the Tevatron experiments, CDF and D0. Combining the measurements of the mass from the 
lepton + jets channel from the two experiments yields a value for the top quark mass which 
has the least error obtained to date. In what follows, we will report on mass measurements in 
three different decay modes of the tt system. The standard model top quark decays into a b 
quark and a W. The final states of the tt system are distinguished by the decay modes of the 
two id’s in the event. These are, 

• the dilepton channel, where both the IH’s decay into a lepton and a neutrino. Due to the 
presence of the two neutrinos, this channel is underconstrained. Both CDF and D0 have 
differing methods to extract the mass from this channel. 

• the all jets channel, where both the IH’s decay into qq pairs. There are no neutrinos in 
this channel that can lead to missing transverse energy. CDF has reported a measurement 
in this channel. 

• the lepton + jets channel, where one W decays into a lepton and neutrino and the other 
decays into qq . Both CDF and D0 have similar methods to extract the mass information. 
The D0 analysis, using multivariate techniques, is new and is being reported for the 
hrst time at this conference. We will describe this channel from both experiments more 
completely than the others, since the global average mass is being determined using this 
channel only and also because the other results have been reported at previous conferences. 

In all that follows, jet energies have been corrected down to the parton level, using Monte Carlo 
models of fragmentation. In making event selection cuts, however, CDF uses uncorrected jet 
energies whereas D0 uses corrected energies. 


2 Kinematic fitting 

In the absence of initial or hnal state gluon radiation, the 6 decay particles of the tt system 
can be described by 18 variables. Each particle is described by a three momentum, its energy 
being determined given its rest mass. The hadronic system that recoils against the tf system is 
described by its pt, which adds another two variables. There are five constraints on the model 
to fit the event, namely that the effective mass of the two W decay particles has to equal the 
W mass, that the effective masses of the top and anti-top decay products have to be equal to 
each other and that the transverse momentum components of the recoiling hadronic system 
have to equal the transverse momentum components of the tt system. This implies that the 
theoretical htting model has 15 free parameters. In the all jets case, all the hnal state particles 
and the transverse momentum components of the recoiling system are observed, yielding 20 
measurements. The system is thus overconstrained by 5, yielding a 5C ht. In the lepton -|- jets 
case, the neutrino three momentum is unknown, yielding 17 measured variables leading to a 2C 
ht. In the dilepton channel, two neutrinos are missing, yielding 14 measured variables. This 
leads to an underconstrained situation (—10 ht). If the top quark mass is specihed, one gets a 



OC case, and the neutrino solutions can be determined for each given top quark mass, leading 
to a likelihood distribution for the top quark mass for each event. Since in general, one does 
not know if a particular jet is the result of a 6 quark decay, there exist several permutations of 
hnal state particles that must be htted for the hypothesis in question. In the all jets case, when 
both the b quarks are tagged, there exist 6 ways of combining the remaining 4 jets into two W 
decay groups. With only one b quark tagged, there exist 10 ways of combining the remaining 
5 jets into the two jets associated with the tagged b jet and another three permutations among 
the remaining three jets to assign the untagged b jet, yielding 30 combinations in all. For 
the lepton + jets channel, there are 12 combinations in the untagged case and 6 combinations 
in the tagged case. For the dilepton channel, there are two combinations for the tagged and 
untagged case. The htting procedure is applied to each combination independently and only 
combinations that meet a goodness of £t criterion are kept in the constrained case. In 
the unconstrained case, the likelihood distributions from the combinations are added up and 
renormalized to obtain the likelihood distribution for the event. When extra jets are present, 
the number of combinations increases rapidly. For this reason, unless otherwise stated, only 
the highest Et jets are used in the £t, the number of jets being the minimum required to fulhll 
the kinematic hypothesis. 


3 The Dilepton channels 

The CDF collaboration p| employs two different techniques in extracting the mass from the 
dilepton channels. Using event selection cuts described in the previous talk |^, and a cut Ht > 
170 GeV, where Ht is the scalar sum of the Et of the various objects, CDF obtains 1 candidate 
in the ee channel, 1 in the /i/i channel and 6 in the e/i channels. The hrst method compares Et 
of the found jets with Monte Carlo from different masses. This relies on the fact the average 
Et of the two b quark jets from top quark decay is directly related to the top quark mass. The 
heavier the top quark, the more energetic are the b quarks on average. This method yields a 
mass of 159 GeV/c^ . 

The second method uses the approximate expression + 2 < > /(I— < 

cos{9ib) >) where is the lepton b quark effective mass and 9ib is the angle between the 
lepton and the b quark. The quantity < cos(6';fc) > is obtained from Monte Carlo and a 
correspondence function between Monte Carlo and data is used to calibrate the result yielding 
a mass of 162 ± 2l.0{stat)1il'^{sys) GeV/c^. 

The D0 collaboration uses a variant of the method proposed by Dalitz, Goldstein and 
Kondo where for each top mass hypothesis, one tries to obtain solutions for the neutrinos. 
This results in 0,2 or 4 solutions. The event four vectors are smeared many times in order to 
estimate the probability of neighbouring events fluctuating to give the observed event. If one 
introduces additional information on the QCD production of the tt system, one can weight 
each solution by a weight proportional to the product of the structure functions and a decay 
probability^. The D0 analysis yields a mass 0 of 158 ± 24.0(sfaf) ± 10(s?/s) GeV/c^. 


4 The All jets channel 

Gnly the CDF collaboration has reported a mass measurement in this channel. CDF performs 
a three constraint ht in this channel, (instead of a possible five constraint ht) by not demanding 
Pt balance between the recoiling system and the tt system. By demanding 5 < Njets < 10 per 
event and the jet TjEt > 200 GeV, CDF observes 142 events with a b quark tag in the silicon 
vertex detector, with a calculated background of 113. This agrees with the rate expected from 


top quark production. Performing a kinematical 3C fit yields a mass |0] of 187± 8 {stat)tl 2 {sys) 
GeV/c^ in this channel. Figure |l| shows the likelihood function for the CDF mass £t in this 
channel. 
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Figure 1; Likelihood function versus mass for the CDF all jets mass determination 


5 The lepton + jets channels 

5.1 CDF results 

The CDF collaboration uses three different techniques to determine the mass from the lepton 
+ jets channels. The hrst is largely unchanged from the time of the top quark discovery |^. 
They require 3 jets with Et >15 GeV and pseudo-rapidity \r]\ < 2.0 and a fourth jet with 
Et >8 GeV and |?7| < 2.4 . They require at least one jet to be 6-tagged with the silicon vertex 
detector or using a secondary lepton. Thirty four events pass the selection criteria. All events 
are treated as a single sample . The estimated non-top background in this sample is 
events. This yields a top quark mass of 175.6 ± 5.7{stat) ± 5.9(sj/s)GeV/c^. 

The second method is called the L** method, where they use in addition the jet charge 
and jet tagging probability to enhance the discriminant. The silicon vertex detector tagging 
probability is used to assign a jet to be of primary (5 quark) or secondary (W decay) origin. 
This information is used to weight each combination. An algorithm is used to estimate the 
leading quark charge to discriminate between b and b jets. This method yields a top quark 
mass of 174.2 ± 5.5{stat) ± 5.3(s|/s)GeV/c^. 

The third method used is termed the optimized method, where they divide the data into 
4 mutually exclusive sub-samples. 

• Events in which a single jet is tagged as a 6 quark jet using the silicon vertex detector 

(SVX). This yields a mass of 176.3 ± 8.2{stat)GeV/c^. 

• Events in which two jets are tagged as b quark jets using the silicon vertex detector (SVX). 

This yields a mass of 174.3 ± 7.9{stat)GeV/. 











• Events in which one or more jets are tagged as a 6 quark using an associated lepton and 
no SVX tag is present. This yields a mass of 140.0 ± 24.1(stot)GeV/c^. 

• Untagged events. This yields a mass of 180.9 ± 6.4(sfaf)GeV/c^. 

The signal to background in the untagged events is increased by demanding all jets to have 
Et > 15 GeV. Background analysis is performed for each individual subsample and separate 
hts are performed for each subsample. The results of the optimized htting method are shown 
in hgure The optimized method results form the main GDF measurement in this channel, 
the other two methods serving as cross checks. The systematic errors for the optimized method 


CDF Preliminary 

M 176.8 ± 4.4 (stat.) ± 4.8 (syst.) GeV/c^ 
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Figure 2: Reconstructed mass for various sub-samples of the GDF optimized lepton -|- jets 
htting. Inset shows the global negative log likelihood function as a function of the top quark 
mass 

have been worked out in detail and are shown in table ^ The sub-sample likelihoods are 
combined into one global likelihood yielding a mass of 176.8 ± 4.4(sfaf) ± 4.8(st/s) GeV/c^ for 
the optimized method. 

5.2 D0 results 

D0 uses two independent multivariate techniques to extract the signal. Multivariate techniques 
permit the separation of signal from background by using an appropriate discriminant that is 
a function of more than one variable. These techniques are superior to the conventional cuts 
method prevalent in high energy physics where signal is separated from background by cutting 
on single variables sequentially. When many variables are needed to separate signal from 
background, the cuts method results in serious losses of signal, especially in cases when the 
amount of signal is small and the signal to background ratio is small in any given variable. 





























GDF Systematic 

Error 

GeV/c2 

D0 systematic 

Error 

GeV/c2 

Soft Gluon + Jet Et scale 

3.6 

Jet energy scale 

4.0 

Monte Garlo event Generators 

1.4 

Generator ISAJET/HERWIG 

3.3 



Generator VEGBGS 

2.6 

Hard Gluon Effects 

2.2 

Multiple interactions 

1.2 

Kinematic and Likelihood fitting methods 

1.5 

LB/NN difference 

1.4 

b-tagging bias 

0.4 

Likelihood fit 

1.0 

Monte Garlo statistics 

0.8 

Monte Garlo statistics 

0.8 

Background spectrum 

0.7 



GDF Total 

4.8 

D0 total 

6.2 


Table 1: Systematic uncertainties for the CDF optimized method and the D0 multivariate 
methods 

5.2.1 Event selection 

D0 selects events with electrons (muons) with Et > 20 GeV, with \ri\ < 2.0(1.7) and missing 
transverse momentum $rp > 20(25) GeV. D0 further demands > 4 jets with Et > 15 GeV 
and \r]jet\ < 2.0. There is also a cut on the scalar Et sum of the W leptonic decay products, 
Et = {Et^^°^ + $t ) > 60 GeV, with \rjw\ < 2.0, for events without a b quark muon tag. Events 
which have a b tag are selected with Pt >4 GeV, with the muon within AR = [Arf + A0^) < 
0.5 of a jet. These cuts yield 91 events of which 7 are b tagged. A 2G top mass fit is performed 
to these events which yields 77 events with < 10, of which 5 are b tagged and ~ 65% are 
background. HERWIG Monte Garlo [§ is used to simulate the signal, VEGBGS MG 0 is used 
to simulate the dominant W + multijet background. The ~ 20% of background events from 
non-W sources are modeled by multijet data, where one of the jets fluctuates to a lepton that 
almost passes the lepton selection. 

5.2.2 Multivariate methods 

D0 dehnes 4 variables to be used in the multivariate analysis, which are so chosen to enable 
us to obtain good signal to background differentiation. These are 

• xi = $T (Missing Et) 

• X 2 = A, the aplanarity, being dehned as 3/2 the least eigenvalue of the normalized 
momentum tensor of the jets and the W boson. 

• Xs = ^ — where Ht = of jets and Hz = E|E^| of the lepton, neutrino and 

the jets, Ez being the momentum component of the object along the beam direction, x^ 
measures the centrality of the event. 

• X 4 ^ = - ^ — , where ARJj'"' is the minimum AR of the six pairs of four jets and E™*"" 

is the smaller jet Et from the minimum AR pair. This variable measures the extent to 
which the jets are clustered together. 

Signal events have larger values of the variables Xi... 0:4 on average than the background events 
that £t with the same value of the top quark mass as the signal. D0 forms two multivariate 















discriminants Dlb and -Dnn using these four variables, where LB stands for the “low bias 
method” and NN denotes a three layer feed-forward neural network with 4 input nodes fed 
by 5 hidden nodes and 1 output node. In the LB method, they hrst parametrize Li{xi) = 
Si{xi)/hi{xi) where Si and bi are the signal and background densities in each variable, integrating 
over the others. A log likelihood In £ = T,iWi In Lj is computed where the weights Wi are adjusted 
slightly away from unity to compensate for any correlation of C with the htted top quark mass. 
Dlb is then dehned as C/{1 -|- C). An event is accepted as a signal event if there exists a b 
quark tag. It is taken as a background event if Dlb <0.43 or if Hrp — < 90 GeV. In the 

neural network approach, the network output Dnn approximates the ratio s(a;)/(s(x) -|- b{x)). 
Figure ^ displays the correlation between Dnn and the htted top quark mass rngt for signal 
events, background events and data. Signal peaks at high values of Dnn at the generated top 
quark mass of 172 GeV/c^ whereas the background peaks at lower values of rngt and Dnn- 
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Figure 3: Events per bin (oc areas of boxes) vs. Dnn (ordinate) and rngt (abscissa) for (a) 
expected 172 GeV/c^ top signal, (b) expected background, and (c) data. 


Figure || shows the distributions of rngt for data (a) passing and (b) failing the LB cut. 
The likelihood analysis proceeds by binning the LB(NN) methods in 40(200) bins in Dlb/nn ^fit 
space and for each bin maximizing the likelihood L{mt,ns,nh) assuming Poissonian statistics 
and Bayes’ theorem|OT; ns,nb are the signal and background events expected in the data. 


Figure ^(c) shows the negative log likelihood thus obtained as a function of the top quark mass. 
This yields a top quark mass of 174.0±5.6(sfaf) GeV/c^ for the LB method and 171.3±6.0(sfaf) 
GeV/c^ for the NN method. Table 0 shows the breakdown of the systematic errors for D0. 
The LB and NN methods are correlated 88 ±4% with each other. The two results are combined 
taking into account these correlations yielding a top quark mass of 173.3 ± 5.6{stat) ± 6.2{sys) 
GeV/c^ . 


5.3 Combining CDF and D0 results 

In order to combine the two results, it is necessary to estimate the common systematic between 
the two experiments. These occur in the areas of Monte Garlo generators for signal and back¬ 
ground, parton fragmentation and luminosity related systematics. While it is certainly not a 
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Figure 4: (a-b) Events per bin vs. mat for events (a) passing or (b) failing the LB cut. 

Histograms are data, filled circles are the predicted mixture of top and background, and open 
triangles are predicted background only, (c) Log likelihood L vs. true top quark mass mt for the 
LB (filled triangles) and NN (open squares) fits, with errors due to finite top MC statistics. 


very well defined process, reasonable people can agree that the common systematic error be¬ 
tween CDF and D0 is conservatively in the neighborhood of 3.0 GeV/c^. With this assumption, 
the CDF optimized method result and the D0 multivariate method result can be combined to 
yield a world average top quark mass of 175.6 ± 5.5 (stat and sys) GeV/c^. Figure ^ shows the 
mass measurements from the various channels and their errors from the two experiments. 


6 Mw vs Mi 


top 


Using a world average W boson mass of 80.410 ± 0.090 GeV/c^[|^, and the currently obtained 


top quark mass, we obtain the comparison shown in hgure Within the currently prevalent 
errors, the standard model is in good agreement, the data perhaps favoring lower Higgs masses. 


7 Conclusion 

We present the measurements of the top quark mass and its world average using lepton -|- 
jets channels from the CDF and D0 experiments at the Fermilab Tevatron. The top quark 
mass measurements using other decay channels are in agreement with the world average within 
errors. A comparison of the top quark mass and W mass world average with the predictions of 
the standard model radiative corrections show no disagreement from what is expected in the 
standard model. Further large improvements in the top quark mass measurement error must 
await data from the upgraded Tevatron and detectors. 
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Figure 5: Summary of mass measurements and errors from the two experiments using various 
methods. The world average, indicated by the vertical line is obtained using the lepton + jets 
channel measurements from CDF and D0 
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Figure 6: Comparison of Mw vs Mtop- Contour shows regions allowed by LEP and SLD data 

























